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Abstract. Recently a stronger statement of Levinson's theorem for the Dirac equa- 
tion was presented, where the limits of the phase shifts at E = itM are related to 



^ ' the numbers of nodes of radial functions at the same energies, respectively. However, 

W . in this letter we show that this statement has to be modified because the limits of 

^ ! the ph ase shifts may be negative for the Dirac equation. 

CN . 

^ ■ Levinson's theorem is one of the fundamental theorem in quantum scattering 

! theory. For the Schrodinger equation, it gives a quantitative relation between the limit 

O I of the phase shift at zero energy 5^(0) and the number of bound states Nf 

5 : Ne = {I/tt} {5,(0} - Seioo)} - {1/2} sm\6e) (1) 

0\ ; 

^ , where the subscript i denotes the angular momentum. As is well know, there is de- 

Q^] generacy of states for the magnetic quantum number due to the spherical symmetry. 

' Usually, this degeneracy is not explicitly expressed in the statement of Levinson's the- 

■ orem. In (1) the phase shifts are determined by comparing them with the phase shift 

^ , . at high energy in order to remove their indetermination of a multiple of 27r. The in- 

determination can also be removed by comparing them with the phase shifts of the 
I free particle, that may be defined as zero^^'^' so that 5^(cxd) in (1) can be removed. In 

^ I fact, the experimental observation for Levinson's theorem shows the jump of the phase 

• shift at low energy in terms of changing the potential t^l. In this letter we will use this 

convention for phase shifts. The second term with the sine square in (1) stands for the 
half bound state, that was first shown by Newton'^1 . 

From the Sturm-Liouville theorem, the number of nodes of the radial function at 
zero energy is equal to the number of the bound states for the Schrodinger equation 
with a non-singular spherically symmetric potential, so that it is related to the phase 
shift through Levinson's theorem. 

Barthelemyt^l first discuss Levinson's theorem for the Dirac equation by the gen- 
eralized Jost function. He stated that Levinson's theorem is valid for positive and 
negative energies separately as in the non-relativistic case. But later this statement 
was found incorrect f^'^l The correct statement of Levinson's theorem for the Dirac 
equation is 

= {I/tt} {6^{M) + 6^{-M)} - {1/2} {sin^(5,(M)) + sm\6^{-M))} (2) 

where k is the standard angular momentum parameter that denotes both the total 
angular momentum j and the orbital angular momentum i: k = j + 1/2 when i = 



1 



J + 1/2, and k. = —{j + 1/2) when i = j — 1/2. Nf^ and ^^(i?) are the number of bound 
states and the phase shift at energy E for the given angular momentum k, respectively. 
The key point of Levinson's theorem for the Dirac equation is that the limits of the 
phase shifts at the thresholds E = ±M (zero momentum) may be positive, zero or 
negative. For example, as the attractive potential becomes strong enough, a scattering 
state of positive energy may change to a bound state, and furthermore, change to a 
scattering state of negative energy. In this case the hmit 6i^{—M) of phase shift at 
E — —M becomes negative. On the other hand, the limit of the phase shift at zero 
energy for the Schrodinger equation only can be zero or positive, if the spherically 
symmetric potential is continuous and finite. 

Recently, Poliatzky'^'^1 transformed the Dirac equation into couple of effective Schrodinger- 
type equations near the thresholds E — ±M and showed a stronger statement of 
Levinson's theorem for the Dirac equation: 

n«(±M) = 5«(±M)/7r - {1/2} sm\6^{±M)) (3) 

where nf^{±M) denote the numbers of nodes for the radial functions with energies ±M, 
respectivelyl^^^ . In Ref. [8] n^{±M) are explained as the numbers of bound states of 
the effective Schrodinger-like equations, that are not easy to count. For the potential 
\V{r)\ < 2M, two numbers are equal to each other due to the Sturm- Liouville theorem. 

A few viewpoints in [8] and [9] have drawn some discussions. Newton'^^' criticized 
the modification of the usual quasi-orthogonality of the scattering wave function. Po- 
liatzky's explanation!^^! that "it is merely a short-hand description of two equations" 
is acceptable. 

Newton's another criticism is that the hmit of phase shift at the threshold cannot 
be 7r/4 in addition to a multiple of tt. In the reply paper^^^l Poliatzky agreed that there 
is a certain constraint to rule out this surprising case'-*^^'. However, he still stated that 
this case cannot be ruled out completely in a conventional proof of Levinson's theorem. 
As a matter of fact, in the proof of Levinson's theorem by the Sturm-Liouville theorem 
the surprising case has been ruled out completely, (see (15b) in Ref. [2] and (14) in 
Ref. [3]) 

It is easy to see that the stronger statement (3) of Levinson's theorem for the Dirac 
equation is not generally correct. The main point is that the limit of the phase shifts 
at thresholds for the Dirac equation may be negative, but the number of nodes, or the 
number of bound states, is a non-negative integer. Why may the limit of the phase shift 
be negative for an effective Schrodinger-type equation? The reason is that the effective 
potential in the effective Schrodinger-type equations is singular when the potential V 
is finite, but strong enough (see Eqs.(22) and (23) in [8] when V = ±2M). For those 
equations Levinson's theorem is not effective. In the same reason, Poliatzky's proof for 
the statement that the sum of the numbers of nodes, n^{M) -|-n«;(— M), is equal to the 
number of bound states is wrong. Besides, there are two radial functions for the 
wave function of the Dirac equation. The numbers of two radial functions may not be 
same. One has to determine which radial function the nodes counted in (3) belong to. 
It seems that Poliatzky assumed'®' that n«;(M) is the number of nodes of mimk, and 
ni^{—M) is that of U2-mk- 

Now, we are going to study this problem in more detail. Discuss the Dirac equation 
with a spherically symmetric potential V{r). For simplicity, we assume the potential 
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is cutoff: V{r) = 0, if r > tq. Let 



KmE 



1 f if^E{r)(t)Km{0,^) 
r \ gKE{r)4>-Km{0, (fi) 



(4) 



where 4>nm{(^-, is the two-component spherical spinor^^l and the radial functions / and 
g satisfy the radial equation: 



I'kE + ~fKE 

9kE ~ ~9kE 



- {E-V + M)g,E, 
{E-V- M)f,E 



(5) 



where the radial functions / and g are proportional to u\ and —U2 used in [8] , respec- 
tively. 

The solution with —k and E can be obtained from that with k and —E by exchange 
f-KE ^ gK-E and by the replacement V — > —V. In the following we only discuss the 
solution with positive k. Assume k > 1 if without notification. Near the origin there 
is only one physically admissible solution so that the ratio of two radial functions at 
r = To is determined. 

In the region r > Tq we have V — There are two oscillatory solutions when 
|£'| > M. Through appropriate combination they can meet the match condition at 
r = ro: 



fnEij-) 



9 Mr) 



= A^{E) 



fnEir) 



9KE{r) 



The phase shift is determined from this condition (see (14) in [3]): 

7r(A;ro)2'^-i A,{E) + 2Mro/(2«; + 1) 



tan5«(£^) ' 
when E > M and E 
tanSKiE) ~ 



{2k - 1)!!(2k - 3)!! A^{E) + 2M{2k - l)/(A;Vo) ' 

M, and 

7r{krof^-^ A^{E) - A;Vo/{2M(2fi; + 1)} 



(6) 



(7a) 



{2k - 1)\\{2k - 3)\\ 



A.{E) 



2k 



2Mro 



{2k - 1){2k - ?>) 



{7b) 



when E < -M and E ~ -M, where k ^ {E'^ - M^Y^'^. Therefore, the phase shift is 
monotonic with respect to the ratio Af^{E) (see (15) in [3]): 



dS,{E) 



dA,{E) 



< when E > M, and 



d5^{E) 



dA,{E) 



> when E < -M (8) 



E 



From (7) we see that the limits of the phase shifts at E = ±M are multiples of tt when 
K > 1. When Af^{M) decreases (or increases) across infinity (5'km('^o) = 0), ^^(M) 
jumps by a TT (or — tt), and when Ai^{—M) decreases (or increases) across the value 
pi, = {2k - l)/(2Mro), 5^{-M) jumps by a -tt (or tt). When k = 1 and Ai{-M) 
decreases (or increases) to the value pi, 5i{—M) jumps by a —n/2 (or tt/2). For this 
critical case where Ai{—M) — pi, there is a finite solution at E = — M, called half 
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bound state (see (16) for k = 1). The case that 5k{M) or 6k{—M) is equal to 7r/4 in 
addition to a multiple of tt is ruled out by (7). 

On the other hand, when |£'| < M there is only one convergent solution in the region 
r > ro so that the match condition (6) is not always satisfied. When the condition 
(6) is satisfied, a bound state appears at this energy. The generahzed Sturm-Liouville 
theorem says (see (17) in [3]) that as the energy E or —V increases, the ratio A^{E) 
decreases monotonically, and as the energy E increases, the ratio fK,E{ro+)/9KE{ro+) 
with < M increases monotonically. Prom the solution in the region r > Tq we have 
(see (18b) in [3]): 



fnEir) 



9KE{r) 



ro+ 

If y = 0, we have (see (18a) in [3]): 

2Mr, 



oo if E = M 

2Mro 



^ = p, .fE = -M (9) 



P2 if E^M 



2k + 1 - (10) 
of E^-M 



It is easy to see from (9) and (10) that there is no overlap for two varying ranges 
of two ratios such that there is no bound state for the case V — 0. When V decreases 
from zero, fKE{ro+) / 9KE{ro+) with \E\ < M docs not change, but A^{E) does decrease. 
Each time Ai^{M) decreases across infinity ((?k_b(^o) = 0), a new overlap for two ranges 
occurs, namely, a scattering state of positive energy becomes a bound state. At the 
same time Si^{M) jumps by tt. On the other hand, each time Ai^{—M) decreases 
across the value pi, an overlap disappears, so that a bound state becomes a scattering 
state of negative energy, 5k{—M) jumps by — tt. If V increases from zero (repulsive 
potential), the process goes conversely, namely, if Ak{—M) increases across the value 
Pi, a scattering state of negative energy changes to a bound state and di^{—M) jumps 
by TT, and if Af^{M) increases across infinity, a bound state changes to a scattering state 
of positive energy and 5^{M) jumps by — tt. For k — 1 we have to consider the possible 
half bound state. 

It is the essence of Levinson's theorem for the Dirac equation. 

Now, we give a calculable example to show the relations between the phase shifts 
and the numbers of nodes explicitly. Consider a square well potential: 

V(r) = I ~ ^ if r < ro . . 

^ ^ \ if r > ro ^ ^ 

The potential V is attractive when A > 0, and is repulsive when A < 0. Solve (5) at 
the energies E = ±M. 
i) E = M 

In the region ro < r where V — 0, the solution for £^ = M is: 



QKMir) = g^M{ro) — 



ro/ ZK+L [yroy Wo 

r 



;i2) 
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In this region there is no node for gnMir), unless gK,M{ro) = 0. In the latter case 
gnuij') = for r > Tq. For this case we say QkM^'^) has a node at r = Tq as usual. 
fKM{r) has one node in the region r > ro if Ak{M) > 0, and no node if Ak{M) < 0. 
Obviously, fnM{f) has a node at ro when Ai^{M) — 0. When A„(M) goes to infinity, 
(gKMiro) = 0) we have a bound state of energy M: 



To 



gnM{r) = 0, for r > ro 



(13) 



Prom the solutions for = M in the region < r < ro we obtain the ratio 
Ai^(M) that decreases monotonically as A increases: 



A.{M) = 



|2M + A 
'2M- |A| 



1/2 



JK+i/2{Piro) 

iJn+l/2{iP2ro) 



J ^«-i/2(«P2ro) ' 
|A|-2M |^/' J.+i/2(p3ro) 

|A| J J«-l/2(p3ro) ' 



if A>0 
if > A > -2M 
if A < -2M 



(14) 



where pi = {|A|(|A| + 2M)}V2, = {|A|(2M - |A|)}V2 and = {|A|(|A| - 2M)}V2, 
and ^^(r) is the Bessel function, 
ii) E^-M 

The solution in the region ro < r for £■ = —M is: 



fn-Mir) = /«_M(ro) 
gK-M{r) = gK-M{ro) 



+ 2^ 



-K+l 



(15) 



In this region there is always a node at infinity for f^^Mir). When fi^^Mijo) = 0, 
fK-uir) = in the region r > rQ. g^^^Mij) has no node if Ak(— M) < pi and one node 
if A^[—M) > pi. When /!«;(— M) = pi we have a bound state of energy —M: 



/«-M(r) = /K-M(ro) — 



gK-M(r) = gK-M{ro) ( — 



-K+l 



(16) 



When K = 1 it is a half bound state. From the solutions for = — M in the region 
r < ro'^^ we obtain the ratio A^{—M) that decreases monotonically as A increases: 



A.{-M) 





1^ 






+ 2M 



A I Y^"^ ^«+i/2(piro) 



A 



1/2 



2M- A 
A-2M 



iJK+i/2{ip2ro) 
JK-i/2{ip2ro) 
4+1/2(^3^0) 
4-1/2 (Psro)' 



if A<0 

if < A < 2M (17) 
if A > 2M 



Carefully counting the nodes of fK±M{f) and gK±M{f) when A increases and de- 
creases from zero, we reach the following conclusions: 
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a) (2) is the correct statement of Levinson's theorem for the Dirac equation. 

b) As a stronger statement of Levinson's theorem for the Dirac equation, (3) is not 
general correct, especially for a very strong potential > 2M). 

c) In the example of a square well potential, (3) can be modified as follows. For 
positive K, if 5^{M) > 0, 5i^{M)/t[ is equal to the number of nodes of gKM{r), which is 
the same as that of fuMir); if ^^(M) < 0, —6i^{M)/tt is equal to the number of nodes 
of g^uir) in the region < r < tq; if 5,{-M) > 0, (5«(-M)/7r) - sm\5,{-M))/2 
is equal to the number of nodes of gn-Aiir); and if M) < 0, — M)/7r) + 
sm'^{6f,{—M))/2 is equal to the number of nodes of g^^Mir) in the region < r < vq 
subtracting the number of its nodes in the region r > tq- However, it is still an open 
problem whether there is and what is the stronger statement of Levinson's theorem for 
the Dirac equation. 

Acknowledgments. The author would like to thank Professor Chen Ning Yang 
for drawing his attention to the papers [8] and [9]. This work was supported by the 
National Natural Science Foundation of China and Grant No. LWTZ-1298 of Chinese 
Academy of Sciences. 

References 

[1] N. Levinson, K. Danske Vidensk. Selsk. Mat-fys. Medd. 25, No.9 (1949). 
[2] Zhong-Qi Ma, J. Math. Phys. 26, 1995 (1985). 
[3] Zhong-Qi Ma, Phys. Rev. D32, 2203 (1985). 

[4] R. G. Newton, Scattering theory of waves and particles, Springer- Verlag, New 
York, 1982, 2nd ed. 

[5] R. G. Newton, J. Math. Phys. 1, 319, (1960). 

[6] M. -C. Barthelemy, Ann. Inst. Henri Poincare, A7, 115 (1967). 

[7] Zhong-Qi Ma and Guang-Jiong Ni, Phys. Rev. D31, 1482 (1985). 

[8] N. Poliatzky, Phys. Rev. Lett. 70, 2507 (1993). 

[9] N. Poliatzky, Helv. Phys. Acta 66, 241 (1993). 

[10] R. G. Newton, Helv. Phys. Acta 67, 20 (1994). 

[11] N. Poliatzky, "Normalization of scattering states and Levinson's theorem: reply 
to a comment by R. G. Newton [Helv. Phys. Acta 67, 20 (1994)]", preprint |hep-| 

I th/941lT97| . 

[12] M. Sassoh de Bianchi, J. Math. Phys. 35, 2719 (1994). 



6 



